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Abstract. The correlation between measured tropospheric 
ozone (O 3 ) and carbon monoxide (CO) has been used ex- 
tensively in tropospheric chemistry studies to explore the 
photochemical characteristics of different regions and to 
evaluate the ability of models to capture these characteris- 
tics. Here, we present the first study that uses multi-year, 
global, vertically resolved, simultaneous and collocated O 3 
and CO satellite (Tropospheric Emission Spectrometer) mea- 
surements, to determine this correlation in the middle/lower 
free troposphere for two different seasons, and to evaluate 
two chemistry-climate models. We find results that are fairly 
robust across different years, altitudes and timescales con- 
sidered, which indicates that the correlation maps presented 
here could be used in future model evaluations. The highest 
positive correlations (around 0 . 8 ) are found in the northern 
Pacific during summer, which is a common feature in the 
observations and the G-PUCCINI model. We make quantita- 
tive comparisons between the models using a single-figure 
metric (C), which we define as the correlation coefficient 
between the modeled and the observed O 3 -CO correlations 
for different regions of the globe. On a global scale, the 
G-PUCCINI model shows a good performance in the sum- 
mer (C = 0.71) and a satisfactory performance in the winter 
(C = 0.52). It captures midlatitude features very well, es- 
pecially in the summer, whereas the performance in regions 
like South America or Central Africa is weaker. The UKCA 
model (C = 0.46/0.15 for July-August/December-January 
on a global scale) performs better in certain regions, such 
as the tropics in winter, and it captures some of the broad 
characteristics of summer extratropical correlations, but it 
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systematically underestimates the O 3 -CO correlations over 
much of the globe. It is noteworthy that the correlations look 
very different in the two models, even though the ozone dis- 
tributions are similar. This demonstrates that this technique 
provides a powerful global constraint for understanding mod- 
eled tropospheric chemical processes. We investigated the 
sources of the correlations by performing a series of sensi- 
tivity experiments. In these, the sign of the correlation is, 
in most cases, insensitive to removing different individual 
emissions, but its magnitude changes downwind of emission 
regions when applying such perturbations. Interestingly, we 
find that the O 3 -CO correlation does not solely reflect the 
strength of O 3 photochemical production, as often assumed 
by earlier studies, but is more complicated and may reflect a 
mixture of different processes such as transport. 


1 Introduction 

Tropospheric ozone (O 3 ) and carbon monoxide (CO) are 
well-known pollutants, majorly affecting the oxidizing ca- 
pacity of the troposphere and exerting a direct (O 3 ) and indi- 
rect (CO and O 3 ) influence on climate. In addition to being 
driven by natural and anthropogenic emissions, their concen- 
trations also depend on a wide range of chemical and physi- 
cal processes that take place in the troposphere (Logan et al., 
1981; Lelieveld and Dentener, 2000). CO is a precursor of 

O 3 , but the relationship between their abundances is not eas- 
ily defined, as it also depends on the concentrations of other 
trace gases and on the meteorological conditions. The strato- 
spheric influence on tropospheric chemistry and the deposi- 
tion of O 3 on the Earth’s surface can make this relationship 
even more complicated. 
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In the past, the correlations between O 3 and CO concen- 
trations from observations have been used in order to under- 
stand the anthropogenic influence on O 3 . The most thor- 
oughly studied area has been the northeastern US and the 
pollution outflow region in the northern Atlantic. There 
have been studies on this region using surface measurements 
(Parrish et al., 1993; Chin et al., 1994; Parrish et al., 1998; 
Cardenas et al., 1998; Li et al., 2002; Honrath et al., 2004; 
Mao and Talbot, 2004) and using aircraft measurements 
(Buhr et al., 1996; Daum et al., 1996; Cooper et al., 2002b,a; 
Huntrieser et al., 2005). They all showed that O 3 -CO cor- 
relations in the summer, when photochemistry is at its peak, 
are positive and often strong, depending on the weather pat- 
terns that prevailed during the measurement period (Cooper 
et al., 2002 a,b). 

Since CO is a good indicator for industrial and biomass 
burning pollution (Logan et al., 1981), a positive correlation 
with O 3 indicates that a region has experienced photochem- 
ical O 3 production from its precursors (including CO). A 
negative O 3 -CO correlation is found in the northern Atlantic 
in the winter (Parrish et al., 1998; Cardenas et al., 1998; Li 
et al., 2002; Mao and Talbot, 2004) which can indicate that 
O 3 destruction through NO x (NO + NO 2 ) titration is favoured 
when pollution increases. However, in Chin et al. (1994) and 
Real et al. (2008) it was also discussed that negative cor- 
relations can sometimes result from secondary CO produc- 
tion from hydrocarbon oxidation (which leads to CO being 
produced but O 3 being consumed), instead of photochemi- 
cal O 3 destruction. Similarly, positive O 3 -CO correlations 
can be explained by CO destruction, rather than O 3 produc- 
tion. However, we expect that these alternative explanations 
do not apply to the majority of cases. We also need to note 
that there are cases when the O 3 -CO correlation is driven by 
factors not affected by chemistry, like when a strong strato- 
spheric intrusion occurs. In such cases, it is expected that the 
correlation will be negative, since stratospheric air is rich in 
O 3 but poor in CO (e.g. Parrish et al., 1998). 

For another pollution outflow region, the northwestern Pa- 
cific, several studies (laffe et al., 1996; Kajii et al., 1998; 
Pochanart et al., 1999; Mauzerall et al., 2000; Tsutsumi and 
Matsueda, 2000; Kato et al., 2004; Suthawaree et al., 2008; 
Tanimoto et al., 2008) found that the broad characteristics 
of the O 3 -CO correlations are similar to those of the north- 
ern Atlantic: positive correlations in the summer and nega- 
tive in the winter, whereas in spring, when stratospheric in- 
trusions are common, correlations are often negative (Hsu 
et al., 2004). Apart from those well-known pollution out- 
flow regions, other correlation analyses using in-situ mea- 
surements have focused on different areas around the globe, 
including continental Asia and America (Chin et al., 1994; 
Naja et al., 2003; Gao et al., 2005; Wang et al., 2006, 2008; 
Ding et al., 2009), Europe (Zahn et al., 2000; Fischer et al., 
2003), the eastern Pacific (Nowak et al., 2004; Bertschi and 
laffe, 2005; Swartzendruber et al., 2008), the tropics (An- 
dreae et al., 1994; Harris et al., 1998; Sanhueza et al., 1999; 


Chan et al., 2002; Takegawa et al., 2003; Wang et al., 2004; 
Henne et al., 2008; Ancellet et al., 2009; Shim et al., 2009) 
and high latitude sites (laffe et al., 1998; Harris et al., 2000). 
These analyses revealed similarities and differences between 
the various regions, which will not be documented here, but 
will be referred to later in this manuscript whenever neces- 
sary. 

Although there have been many studies that examined 
the important topic of O 3 -CO correlations, their reliance on 
surface/aircraft observations has left much of the globe un- 
explored. The aircraft campaigns provided relatively good 
(though far from complete) spatial coverage, but for limited 
periods. The satellite observations from the Tropospheric 
Emissions Spectrometer (TES; Beer et al., 2001) provide a 
unique opportunity to study the global multi-year O 3 -CO 
correlation from collocated and simultaneous, vertically re- 
solved tropospheric O 3 and CO data. There have been a few 
recent studies that have taken advantage of this dataset for 
this purpose, but they have either focused on particular re- 
gions (Zhang et al., 2008; Shim et al., 2009) or examined a 
limited time-period (Zhang et al., 2006). Here, we expand 
the scope of Zhang et al. (2006) (who focused on July 2005) 
to examine the middle/lower free tropospheric O 3 -CO cor- 
relation patterns around the globe in two different seasons 
(Northern Hemisphere summer/winter), using a multi-year 
(2005-08) dataset. Since the O 3 -CO correlation analysis pro- 
vides information on whether a model captures O 3 levels for 
the right reasons (Chin et al., 1994), we also employ two 
different models that we evaluate in terms of their ability to 
capture the correlation patterns. Furthermore, we perform a 
sensitivity analysis, in which the effect of different emissions 
on the correlations is examined. 

In Sect. 2 we provide some information about the TES ob- 
servations and the two models used (G-PUCCINI, UKCA). 
Section 3 presents CO, O 3 and O 3 -CO correlations as ob- 
served by TES and as simulated in the models, as well as an 
analysis of performance of the models by using the correla- 
tion coefficient between the modeled and the observed O 3 - 
CO correlations (C). Section 4 analyzes the results of a sen- 
sitivity analysis that examines whether emissions drive the 
broad features of the correlations and Sect. 5 discusses poten- 
tial reasons for differences between the models and TES ob- 
servations. Finally, the conclusions are presented in Sect. 6 . 
The purpose of this study is not to separately analyze and 
explain the detailed features of correlations in all different 
parts of the globe, but to document global correlation maps 
that can be used for future model evaluation, and discuss the 
broad features found, with an emphasis on the possible role 
of emissions as a driver. 
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2 Datasets 

2.1 TES satellite data 

The Tropospheric Emission Spectrometer (TES) is a high- 
resolution (0.1 cm -1 ), infrared, Fourier Transform spectrom- 
eter aboard the NASA Aura satellite. Aura satellite follows 
a polar Sun-synchronous orbit with an equator crossing time 
at 01:45 and 13:45 local time and has a repeat cycle of 16 
days. TES covers a wide spectral range (650 to 3050 cm -1 ) 
and has an averaged nadir footprint of about 5 km by 8 km 
(Beer et al., 2001). The spectral radiances measured by TES 
are used to retrieve the species’s (e.g. O3, CO) mixing ratio 
profiles using the optimal estimation method (Rodgers, 2000; 
Bowman et al., 2002, 2006). The retrieved natural logarithm 
of the mixing ratio profile x can be expressed as: 

x = x a + A(x-x a ) + e (1) 

where x and x a are the natural logarithms of the true atmo- 
spheric profile and of the a priori profile respectively, A is the 
averaging kernel matrix, expressing the sensitivity of the re- 
trieval to the true state, and e is the observational error, x, x a , 
A and e are part of the standard TES product. Also, a master 
quality flag that we use to exclude bad data points is avail- 
able with the TES product. x a consists of monthly means 
produced from MOZART model simulations. The TES prod- 
ucts are on 67 vertical levels with a varying layer thickness. 
In this study we use the version 4 of TES global survey data, 
each of which consists of data from 16 full orbits. 

TES O3 data have been evaluated by comparison to 
ozonesondes (e.g. Worden et al., 2007; Nassar et al., 2008; 
Bowman et al., 2009; Boxe et al., 2010), aircraft data (e.g. 
Richards et al., 2008), and satellite data (OMEMLS) (e.g. 
Osterman et al., 2008). These studies show that TES O3 
generally has a positive bias that varies between 3-10 ppbv 
(Nassar et al., 2008) in the troposphere. TES CO has been 
validated against aircraft data and has been shown to have a 
slightly negative (<10%) bias in midlatitudes and a slightly 
positive bias (<10%) in the tropics (Luo et al., 2007; Lopez 
et al., 2008). The TES O3 and CO data have been used ex- 
tensively to examine a variety of scientific topics (e.g. Par- 
rington et al., 2008; Bowman et al., 2009; Shim et al., 2009; 
Hegarty et al., 2010), including the investigation of the 3-D 

03 radiative forcing as observed from space (Aghedo et al., 

2011 ). 

2.2 G-PUCCINI model data 

G-PUCCINI is the Goddard Institute for Space Studies 
(GISS) model for Physical Understanding of Composition- 
Climate INteractions and Impacts (Shindell et al., 2006b). 
The composition model (including both gases and aerosols) 
is fully embedded in the GISS modelE climate model 
(Schmidt et al., 2006). Tropospheric and stratospheric chem- 
istry are simulated simultaneously. Tropospheric chemistry 


includes basic NO x -HO x -O x -CO-CH4 interactions as well as 
PAN, isoprene, alkyl nitrates, aldehydes, alkenes and paraf- 
fins. The lumped hydrocarbon family scheme was derived 
from the Carbon Bond Mechanism-4 (Gery et al., 1989) and 
from the more extensive Regional Atmospheric Chemistry 
Model, following Houweling et al. (1998). The stratospheric 
chemistry includes chlorine- and bromine-containing com- 
pounds, and CFC and N2O source gases. The full scheme 
includes 155 chemical reactions among 55 species. The 
main additions to the chemistry of the previous versions of 
G-PUCCINI are (a) the addition of acetone to the hydro- 
carbons following Houweling et al. (1998), (b) the depen- 
dence of polar stratospheric cloud formation on the abun- 
dance of nitric acid, water vapor and temperature (Hanson 
and Mauersberger, 1988), and (c) the addition of a reac- 
tion pathway for HO2+NO to yield HNO3 (Butkovskaya 
et al., 2007). Photolysis rates are calculated using the Fast-J2 
scheme (Wild et al., 2000; Bian and Prather, 2002), whereas 
other chemical reaction rate coefficients are from JPL-2000 
(Sander et al., 2000). Tracer transport uses a non-diffusive 
quadratic upstream scheme (Prather, 1986). A full descrip- 
tion of the model is given in Shindell et al. (2006b) and ref- 
erences therein. 

The G-PUCCINI model has been extensively compared 
with observations and other models (e.g. Stevenson et al., 
2006; Dentener et al., 2006; Shindell et al., 2006b,a; Sander- 
son et al., 2008; Fiore et al., 2009; Lamarque et al., 2010). 
It performs well at simulating tropospheric O3 and CO con- 
centrations, with some strengths and weaknesses that will be 
highlighted later in the discussion. 

We ran the model at a 2° latitude by 2.5° longitude hori- 
zontal resolution, with increased effective resolution for trac- 
ers by carrying higher order moments at each grid box. This 
configuration has 40 vertical hybrid pressure layers from the 
surface to 0.1 hPa (?»80km). Simulations were performed 
using observed sea-surface temperatures (Rayner et al., 
2003) and linear relaxation of winds toward NCEP/NCAR 
reanalysis (Kalnay et al., 1996). We use present-day an- 
thropogenic emissions for a variety of sectors produced in 
support of the Intergovernmental Panel on Climate Change 
(IPCC) Fifth Assessment Report (AR5) (Lamarque et al., 
2010). For biomass burning we use data from the Global 
Fire Emissions Database version 2 (GFED2) (Randerson 
et al., 2007). Isoprene emissions are a function of vegeta- 
tion type and leaf area index with responses to temperature 
and solar radiation based on the algorithm of Guenther et al. 
(1995, 2006). Vegetation emissions of alkenes and paraffin 
come from the GEIA dataset, and are based on Guenther et 
al. (1995). Lightning NO x emissions depend on the climate 
model’s convection based on the parameterization of Price 
et al. (1997). The abundances of long-lived species (includ- 
ing methane) have been prescribed to observations. 
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2.3 UKCA model data 

UKCA is the UK Chemistry Aerosol Model (Morgenstern 
et al., 2009; O’Connor et al., 2011), which is based upon 
the UK Met Office’s Unified Model (UM) (Staniforth et al., 
2005). We employ the model set up of Telford et al. (2010) 
who used the tropospheric version of the model (O’Connor 
et al., 2011). This has a medium-sized chemistry, simulat- 
ing the NO x , HO x and O x chemical cycles and the oxida- 
tion of CO, ethane, propane and isoprene, with 132 reactions 
among 60 species. Long-lived tracers, such as methane, are 
set to constant values throughout the atmosphere and O3 and 
NO x are overwritten above 30hPa with zonal mean values 
from the Cambridge 2-D model (Law and Pyle, 1993). Pho- 
tolysis rates are calculated using the offline scheme of Law 
et al. (1998). Tracer transport is performed using the semi- 
Lagrangian scheme of Priestley (1993). An evaluation of the 
tropospheric version of the model can be found in O’Connor 
et al. (2011). 

We ran the model at a 2.5° latitude by 3.75° longitude 
horizontal resolution. This configuration has 60 vertical hy- 
brid height layers from the surface to 84 km. Simulations 
were constrained by sea-surface temperatures (Rayner et al., 
2003) and with winds and temperature “nudged” (Telford 
et al., 2008) toward ECMWF Operational analyses. We 
use the IPCC AR5 emissions inventory from 2000 for both 
the anthropogenic and biomass burning emissions (Lamar- 
que et al., 2010). Biogenic emissions of isoprene are from 
Guenther et al. (1995) and there are 40 Tg of acetone emit- 
ted per year, spatially distributed as the isoprene emissions 
(O’Connor et al., 2011). In addition, there are 45 Tg of CO 
emitted from the oceans and 18.4 Tg of NO emitted from 
soils with distributions taken from the GEIA project 1 . Light- 
ning NO x emissions depend on the climate model’s convec- 
tion based on the parameterization of Price and Rind (1994). 

2.4 Data post-processing 

Similarly to when comparing TES data with in-situ measure- 
ments (sondes, aircraft; e.g. see Jones et al. (2003); Worden 
et al. (2007); Shim et al. (2009)), we have post-processed 
the model output in order to be directly comparable with the 
TES observations. This required: (i) sampling of the model 
3-hourly output at the times and locations of the TES mea- 
surements, (ii) interpolating model data onto the 67 TES ver- 
tical pressure levels, and (iii) application of the TES a pri- 
ori profiles and averaging kernels (hereafter, “TES opera- 
tors”), in order to account for the measurement sensitivity 
(see Sect. 2.1) at different times and locations. Step (iii) in- 
volves the use of Eq. (1), but with the logarithm of the mixing 
ratios from the model (extracted and vertically interpolated) 
now replacing the true atmospheric profile. When the aver- 
aging kernel goes to zero, the model post-processed output 

1 http://www.geiacenter.org/inventories/present.html 


will revert to the a priori profile values, as would the TES 
measurements. TES has been shown to be relatively sensi- 
tive in the vertical region on which we are focusing in this 
study (800-400 hPa), in several different geographical areas 
(Zhang et al., 2006; Worden et al., 2007; Richards et al., 
2008; Osterman et al., 2008; Bowman et al., 2009; Shim 
et al., 2009; Hegarty et al., 2009). 

Note that we have binned the processed observational and 
model data onto a 4° latitude by 5° longitude grid in order 
to smooth-out gaps in the observations. A similar smoothing 
was done in the correlation analysis of Zhang et al. (2006), 
but on a coarser grid (10° x 10°), since they considered only 
one month. 

3 Results from TES and models 

Figures 1 and 2 show the mean concentrations of CO and 
O3 for the period of interest (2005-08) and for the two sea- 
sons of interest (July-August and December-January, re- 
spectively). From Fig. 1, it is obvious that both models are 
able to capture the CO spatial variability fairly well: the 
northern hemispheric concentrations are much higher, and 
the biomass burning peaks are clearly visible for both sea- 
sons. Also, CO is more abundant in the Northern Hemi- 
sphere winter, which is reasonable due to the lower amounts 
of hydroxyl radical (OH) that are present in the troposphere 
in that season (reduced oxidizing capacity). The average CO 
concentrations in both models are higher than in TES for both 
seasons. This is in contrast to earlier results from a multi- 
model study (Shindell et al., 2006a), which found that the 
majority of models under-predicted CO in the troposphere, 
especially at boreal latitudes. Our result suggests that proba- 
bly the new anthropogenic emissions dataset that we are us- 
ing (AR5) includes improved emission estimates. Also, the 
differences between modeled and measured CO may have to 
do with the low biases of TES CO (see Sect. 2.1) especially 
at midlatitudes. 

The spatial distribution of O3 in TES is captured well by 
the models (Fig. 2): there is high O3 in anthropogenically 
polluted and biomass burning regions, and over the former, 
O3 is highest in the summer, when photochemical activity is 
at its peak. Also, the highest values of ozone in the sum- 
mertime are seen over Middle East, North Africa and central 
Asia in all datasets, consistent with previous studies for the 
middle troposphere (Liu et al., 2009; Worden et al., 2009). 
The G-PUCCINI model captures the levels of O3 very well, 
in agreement with what Shindell et al. (2006b) found for the 
middle troposphere. In the UKCA model, the spatial pattern 
is captured well, though the O3 levels appear to be somewhat 
lower than the observed. However, note that TES O3 has 
been found to be positively biased (see Sect. 2.1) in much of 
the globe, which implies that the UKCA is likely doing well, 
and better than this comparison implies. 
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Fig. 1 . Mean 2005-08 middle/lower free tropospheric CO concentrations for TES (a, b), G-PUCCINI (c, d) and UKCA (e, f) in July-August 
(left panels) and December-January (right panels). We used data from 7 TES pressure levels between 800 and 400 hPa. TES sampling and 
operators have been applied to the model output. The data have been smoothed by averaging on a 4x5 degree grid. 


Figure 3 shows global maps of the correlation coefficients 
between the 2005-08 O 3 and CO time series (daily means), 
as observed from TES and as simulated by the two models, 
for the two seasons of interest. In the TES map for July- 
August (panel a), it can be seen that several regions have 
positive correlations, with the highest ones being in the west- 
ern part of the northern Pacific (up to 0.8). Positive sum- 
mer O 3 -CO correlations were also found in in-situ measure- 
ments within this region (e.g. Kajii et al., 1998; Mauzerall 
et al., 2000; Suthawaree et al., 2008). In our study, posi- 
tive but lower correlations also extend further East, towards 
the US. In fact, the positive correlations are also found over 
the US, and become reinforced over the other major midlati- 
tude pollution outflow region, from the east US coast towards 
the northern Atlantic, where we find positive correlations of 


up to 0.4-0. 6 . Positive summer O 3 -CO correlations in this 
area are also in agreement with previous studies (e.g. Parrish 
et al., 1993, 1998; Honrath et al., 2004). 

Other regions that stand out as areas of positive O 3 -CO 
correlations include a big part of the tropics and subtrop- 
ics, especially close to biomass burning regions. These ar- 
eas have not been examined in detail in terms of the O 3 -CO 
correlation in the past, but two studies for spring by Collins 
et al. (1996) (southern Atlantic) and by Shim et al. (2009) 
(Gulf of Mexico) revealed positive correlations. Also, pos- 
itive correlations were found during fire events in northern 
Australia (Takegawa et al., 2003) and at a mountain site in 
Kenya (Henne et al., 2008), though in the latter correlations 
were weak. 
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Fig. 2. Mean 2005-08 middle/lower free tropospheric O 3 concentrations for TES (a, b), G-PUCCINI (c, d) and UKCA (e, f) in July-August 
(left panels) and December-January (right panels). We used data from 7 TES pressure levels between 800 and 400 hPa. TES sampling and 
operators kernels have been applied to the model output. The data have been smoothed by averaging on a 4x5 degree grid. 


Areas of negative correlations in July-August include high 
northern latitudes and the deserts of the Sahara and the Mid- 
dle East (regions not studied in terms of O 3 -CO correlation 
in the past). Also, much of continental Asia shows negative 
or weak correlations, in agreement with Naja et al. (2003); 
Wang et al. (2006), while areas closer to the large urban cen- 
ters of East Asia show mainly positive correlations, in agree- 
ment with Gao et al. (2005); Wang et al. (2008); Ding et al. 
(2009). Apart from qualitatively agreeing fairly well with 
past in-situ observational studies, the TES O 3 -CO correlation 
patterns that we find also broadly agree with the findings of 
Zhang et al. (2006) for July 2005 (where they use TES and 
the GEOS-Chem model). 


The positive O 3 -CO correlations found in TES data for 
the midlatitude pollution outflow regions in the summer, are 
still found in the winter (Fig. 3b), but with less latitudinal 
coverage, extending mainly up to 30-40° N. The higher lat- 
itude areas in northern Atlantic/Pacific now show negative 
correlations, in agreement with previous studies (e.g. Par- 
rish et al., 1998; Kajii et al., 1998; Li et al., 2002; Suthawa- 
ree et al., 2008). This indicates that these latitudes, due to 
the lack of sunlight, do not experience intense photochem- 
istry and are dominated by O 3 destruction and/or dry deposi- 
tion. Occasionally, stratosphere-troposphere exchange (STE) 
can also drive these negative correlations, which, however, is 
more frequent in spring (Zhang et al., 2008). The tropics 
show positive correlations in regions where wildfire burning 
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Fig. 3. Correlation coefficient of daily mean 2005-08 middle/lower free tropospheric O 3 and CO concentrations for TES (a, b), G-PUCCINI 
(c, d) and UKCA (e, f) in July-August (left panels) and December-January (right panels). We used data from 7 TES pressure levels between 
800 and 400 hPa. TES sampling and operators have been applied to the model output. Also, the observational error has been taken into 
account for the model correlations, following Eq. (3) of Zhang et al. (2006). Stippling indicates the areas for which less than 40 sample 
members (i.e. days) were available for the correlation calculation. These are the same (but not shown) in the middle and lower panels (c, d, 
e, f), since the geographical and temporal sampling was identical with the upper panels. Within these regions with >40 sample members, 
correlations are statistically significant only where they are higher than 0.3 or lower than —0.3, based on a Pearsons test (95 % confidence 
interval). The data have been smoothed by averaging on a 4x5 degree grid. 


is expected during the December-January season, like West 
Africa and a large part of Central Africa. Most of the re- 
gions that showed negative correlations in July-August, still 
have negative correlations in December-January, like, for in- 
stance, the northern high latitudes, in agreement with Harris 
et al. (2000) who studied O 3 -CO correlations in Alaska dur- 
ing winter. 

Apart from the O 3 -CO correlation coefficients, we 
also calculated the slopes of the TES O 3 -CO regres- 
sion line, which represent the O 3 -CO enhancement ratio 
(c/[ 03 ]/< 7 [C 0 ]) and are shown in Fig. 4a, b. Slopes are useful 


to indicate the strength of the O 3 -CO relationships, and as 
such they have been reported in several observational stud- 
ies (e.g. Parrish et al., 1998; Li et al., 2002). It is clear that 
the patterns of the slopes follow those of the correlations 
(Fig. 3a, b) very well. Additionally, in order to understand 
which regions with high correlations are associated also with 
high absolute variability of pollution, we calculated the TES 
O 3 -CO covariance (see Fig. 4c, d), from the same data used 
for the calculation of the correlation. It is evident that the ar- 
eas with highest covariance are the biomass burning regions 
in both seasons and the midlatitude pollution outflow regions 
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Fig. 4. O 3 -CO slope (c?[03]/rf[CO] ) (a, b) and covariance (c, d) from the 2005-08 800-400 hPa TES data for July-August (a, c) and 
December-January (b, d). 


in the Northern-Hemisphere summer. On the other hand, the 
regions with the lowest absolute covariances cover most of 
the southern middle/high latitudes (>30° S) and much of the 
northern high latitudes (>50° N), suggesting that any dis- 
crepances that we have seen so far regarding the correlation 
in these regions are not important. 

For the calculation of the model O 3 -CO correlations (see 
Fig. 3c, d, e, f), we have applied the observational error pro- 
vided from the TES datasets to the model output using the 
following equation (Zhang et al., 2006): 


Corr (xjXj) 



Cov(.rV.) 

1 J 


+ x?[S] iiy j* 2 x ,+x' j 2 [S]jj 


( 2 ) 


where Corr (xtxj) is the correlation between any two ele- 
ments x i and Xj of the retrieved model O 3 and CO profiles, 
x) and x'f are the same elements under noise-free conditions, 

Cov(xV ; ) is the covariance of these elements, x' ; and x\ are 
1 J 1 J 

their mean values, oy and ay are their standard deviations, 

and finally S is the spectral measurement error covariance 

(we assumed that spectral measurement errors are small in 

the conversion from log mixing ratio to mixing ratio, and that 

the spectral measurement error covariance is the same for all 

scenes, as in Zhang et al., 2006). From the above formula it is 


clear that the correlations in the raw model data are expected 
to be higher than those shown in Fig. 3, but the application 
of the observational error reduced them and made the results 
directly comparable with TES. This was something that we 
did not need to do for Fig. 1 and Fig. 2, since on a monthly 
basis the error is expected to be zero. Note that we performed 
tests in which the observational error was varied by ±25 %, 
to examine whether the model results depend on this quan- 
tity, but found that the main features are fairly insensitive to 
those perturbations. 

G-PUCCINI model results capture many of the correlation 
patterns seen in the observations remarkably well for both 
seasons (Fig. 3c, d). There are certain features that are not 
captured, such as when the correlations become negative in 
the winter at high-latitudes over the Atlantic, but generally, 
many of the characteristics are simulated, including the sea- 
sonal variation of correlations over the northern Pacific. The 
UKCA model (Fig. 3e, f) map shows some similarities with 
TES in terms of distribution of positive and negative corre- 
lations. The main difference is that the correlations are sys- 
tematically lower than the observations for both seasons. 

In Fig. 5, we examine the relationship between the mod- 
eled and observed O 3 -CO correlations using data from all 
over the globe. It is evident that the G-PUCCINI model does 
well in capturing the summer O 3 -CO correlations, while in 
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Fig. 5. Scatter plots of modeled versus observed O 3 -CO correlations for July-August (a, c) and December-January (b, d). The data that 
have been used are the same as for Fig. 3, and have undergone the same processing. 


the winter there seems to be a slight tendency for overesti- 
mation. The UKCA model tends to underestimate the O 3 - 
CO correlations in both seasons, especially in the summer. 
Additionally, in order to quantitatively compare the models, 
we calculate the correlation coefficient (C) between modeled 
and observed O 3 -CO correlations for different regions of the 
globe in both seasons (Table 1). It is clear that the areas 
where the G-PUCCINI model does particularly well cover 
much of the northern midlatitudes in the summer, with C 
higher than 0.6 for North America, northern Atlantic, Eu- 
rope and west northern Pacific. Areas with less good per- 
formance can be found in the tropics (e.g. South America, 
Central Africa), especially in the winter. The performance 
of the UKCA model is satisfactory for the northern midlat- 
itudes in the summer, e.g. exhibiting a good performance at 
the Pacific outflow region. However, with the exception of 
the east northern Pacific, the performance at midlatitudes is 
not good in the winter. On the other hand, there are specific 
cases where the UKCA model seems to have more skill on 
capturing the O 3 -CO correlation, such as in remote regions 
like South America and the east tropical Pacific in the winter. 


4 Sensitivity analysis 
4.1 Robustness of results 

As discussed in Sect. 3, to calculate the model correlation 
coefficients shown in Fig. 3 (c, d, e, f), we have taken into 
consideration the observational error, in order to have a more 
consistent comparison with TES correlations. In Fig. 6 a, we 
show how the G-PUCCINI model July-August correlations 
would look, had we not accounted for this error. As expected, 
the sign of correlations is identical with Fig. 3c, but the mag- 
nitude of correlations is larger. We underline that the magni- 
tude of correlations seen in Fig. 6 a is expected to be closer 
to what could be found in reality in the troposphere, since the 
correlations have not been degraded by observational error. 

In Fig. 6 b, we show how the G-PUCCINI O 3 -CO corre- 
lations would look if the TES operators had not been ap- 
plied to the model output (so that the only post-processing 
would have been the TES time/location sampling). It is clear 
that the results are very similar to Fig. 6 a, with the excep- 
tion of the Southern Ocean, where correlations change from 
negative to positive. This indicates that this region was heav- 
ily weighted towards the TES a priori. However, as seen in 
Fig. 4c, d, this is an area with very low O 3 -CO covariance, 
thus this difference in correlations when ignoring the TES 
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Table 1 . Performance (using our C metric) in terms of capturing the O 3 -CO correlations, of the G-PUCCINI and UKCA models at various 
regions of the globe for two seasons (July-August and December-January). C represents the correlation coefficient between the modeled 
and the TES O 3 -CO correlations. The figures in italics are not statistical significant based on a Pearsons test (95 % confidence interval). Note 
that we have avoided using regions that generally have low data availability (e.g. North Africa) based on Fig. 3a. b. 


G-PUCCINI UKCA G-PUCCINI UKCA 
(Jul+Aug) (Jul+Aug) (Dec+Jan) (Dec+Jan) 


Global 

0.71 

0.46 

0.52 

0.15 

North America (120° W-70° W; 35° N-65° N) 

0.70 

0.33 

0.32 

-0.49 

North Atlantic (70° W-15° W; 35° N-65° N) 

0.81 

0.53 

0.57 

-0.39 

Europe (10° W-35° E; 35° N-65° N) 

0.67 

0.43 

0.31 

- 0.20 

Siberia (60° E-130° E; 45° N-70° N) 

0.15 

0.28 

0.10 

0.16 

West Northern Pacific (130° E-180° W; 20° N-50° N) 

0.72 

0.64 

0.58 

0.24 

East Northern Pacific (180° W-120° W; 20° N-50° N) 

0.40 

0.30 

0.68 

0.50 

South America (80° W-35° W; 30° S-10° N) 

0.23 

0.27 

0.18 

0.40 

South Atlantic (30° W-0° E; 30° S-0° N) 

0.35 

0.04 

0.22 

- 0.04 

Central Africa (10° E-40° E; 20° S-5° N) 

0.17 

0.07 

- 0.04 

0.23 

Indian Ocean (50° E-90° E; 30° S-10° N) 

0.62 

- 0.07 

-0.31 

- 0.15 

West Tropical Pacific (100° E-150° E; 20° S-20° N) 

0.40 

0.24 

0.17 

0.28 

East Tropical Pacific (130° W-85° W; 20° S-20° N) 

0.41 

0.06 

0.36 

0.58 


operators is not important, as it is not associated with large 
enough O3 and CO fluctuations. Figure 6c shows how the 
correlations would look, had there been no TES time/location 
sampling. The plot indicates that the features of the maps that 
we have shown so far do not depend on the TES sampling. 

Figure 6d shows the O3-CO correlations for a single pres- 
sure level (562 hPa) within the vertical region of interest 
(800-400 hPa), rather than for the 80(.M-00 hPa mean. One 
can conclude that our results do not depend on the choice of 
height in the middle/lower free troposphere. The results ap- 
pear to change more when looking at the upper troposphere 
(see panel (e) for correlations at 316 hPa), with e.g. correla- 
tions getting significantly lower over northern high latitudes, 
most likely due to the influence of STE. However, there are 
several features, especially in the northern midlatitudes, that 
remain similar even at this altitude. 

From Fig. 6f, it can be seen that our results do not only 
apply for day-to-day variations of O3 and CO, but are sim- 
ilar when considering somewhat longer timescales (quasi- 
weekly). Note that we have also tested how the correlations 
change when taking into account the raw 3-hourly model 
data instead of the daily averages, and how they change when 
deseasonalizing the O3/CO time series (not shown). Again, 
little change has been found in our results (expect for de- 
seasonalization effects in the Southern Ocean, which were 
more drastic, but only when calculating correlations using 
days from all months of the year). Later in the manuscript 
(Fig. 8a and Fig. 9a), we also discuss that the O3-CO corre- 
lations are not sensitive to the year that we choose. 

For the UKCA model (Fig. 7), similar conclusions about 
the robustness of the results can be drawn. There are some 
notable exceptions, like the fact that when not applying the 
TES operators to the UKCA post-processing (Fig. 7b), (i) 


the northern high-latitude (>60° N) correlations become al- 
most entirely negative, which is more in agreement with the 
observations, and (ii) the positive correlations over biomass 
burning regions in Africa become more widespread, which 
may be related to the fact that the convective processes are 
different in the two models. This suggests that the observed 
correlations can provide a useful constraint for vertical mix- 
ing, something that can be examined further in future work. 

However, in general the correlation features found in our 
analysis seem to be robust and to characterize each model. 
Note that the same robustness exists in the observational 
data (no dependence on year, altitude within 800-400 hPa, 
timescale etc.). Additionally, most of the similarities found 
across different estimates, as described in Fig. 6 and Fig. 7 
for July-August, are also present in G-PUCCINI/UKCA re- 
sults for December-January. In Fig. 6 g and Fig. 7 g we 
show the 2005-08 December-January correlations from the 
raw model output (compare with Fig. 3d and Fig. 3D just for 
reference. Despite the similarities of the results from all the 
different ways of estimating the O3-CO correlations within 
each model, the results from the two different models remain 
significantly different. 

4.2 Perturbing emissions 

Next, we analyze the results of four sensitivity model runs 
that were performed in order to assess the role of emissions 
in driving the correlation patterns. We focus on July-August, 
as we have demonstrated that this is the period with the most 
intense features in the Northern Hemisphere. Through tests 
that we have performed by comparing the O3-CO correla- 
tions of 2005-08 with those from each individual year (e.g. 
compare Fig. 8a with Fig. 6c) we have found that the choice 
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Fig. 6. Testing the robustness of G-PUCCINI 2005-08 correlation coefficients seen in Fig. 3: (a) same as Fig. 3c but without the observational 
error; (b) same as (a) but without the TES operators application; (c) same as (b) but without the TES sampling (i.e. raw model output); (d) 
same as (c) but for 562 hPa instead of 800-400 hPa average; (e) same as (c) but for 316hPa instead of 800-400 hPa average; (f) same as 
(c) but calculated using 5-daily (quasi- weekly) averages instead of daily. Also, we show the correlations from the raw model output for 
December-January in (g) (compare with Fig. 3d). 


of year does not matter much. So we performed the sensitiv- 
ity experiments only for year 2006. 

In each of the runs, one type of emission is “switched 
off” (no_ANTH for no anthropogenic, no_BB for no biomass 
burning, noXIGHT for no lightning and no_VEG for no veg- 
etation emissions). Figure 8a and Fig. 9a show the results 
from the unperturbed runs (BASE) in 2006. 

In both models, it is clear that biomass burning is the emis- 
sion that has the highest impact on the O3-CO correlation in 
the tropics (Fig. 8c, Fig. 9c). Downwind of Central Africa 


and of South America, there is even a change in the sign of 
the correlation coefficients (from positive to negative). This 
implies that the photochemistry of these regions heavily de- 
pends on biomass burning. We expect that in December- 
January, the large impact of biomass burning will be seen 
in north equatorial Africa, where biomass burning is intense 
during that season. It is interesting to note that biogenic 
emissions (which are very important in the tropics) are not 
enough to sustain positive correlations in a large part of these 
regions. There appear to be some regions of positive changes 
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Fig. 7. Testing the robustness of UKCA 2005-08 correlation coefficients seen in Fig. 3: (a) same as Fig. 3e but without the observational 
error; (b) same as (a) but without the TES operators application; (c) same as (b) but without the TES sampling (i.e. raw model output); 
(d) same as (c) but for 562 hPa instead of 80(M-00hPa average; (e) same as (c) but for 316hPa instead of 800-400 hPa average; (f) same 
as (c) but calculated using 5-daily (quasi-weekly) averages instead of daily. Also, we show the correlations from the raw model output for 
December-January in (g) (compare with Fig. 3f). 


in the correlation coefficient in both models, like the outflow 
regions of pollution from Indonesia (Indian Ocean and equa- 
torial Pacific). These are mostly areas with positive and in 
some cases high correlations, which become even higher fol- 
lowing this perturbation. 

The effect of the anthropogenic emissions (Fig. 8 b, 
Fig. 9b) is more straightforward, especially at midlatitudes, 
where we find mostly reductions in the O 3 -CO correlations 
when removing them. An exception, in the G-PUCCINI 
model, is some parts of East Asia, which possibly are dom- 


inated by NO x titration, where correlations are increased in 
the no_ANTH run. However, this feature is not seen in the 
UKCA model. Note that in the 110 A NTH run, correlations 
preserve their sign in most of the regions, indicating that 
there are still other types of emissions in key areas (e.g. bio- 
genic, for the north Atlantic) that preserve the O 3 production 
processes. In the no_VEG run (panels d), correlations are 
again reduced in regions like, for example, the north Atlantic 
Ocean, but do not change sign. There are regions where pos- 
itive correlations increase, like in the west coast of South 
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Fig. 8. Correlation coefficient of daily mean July-August 2006 middle/lower free tropospheric O 3 and CO concentrations in the BASE 
G-PUCCINI run (unperturbed; a) and in 4 sensitivity runs in summer. The four runs are: (b) no anthropogenic emissions, (c) no biomass 
burning emissions, (d) no emissions from vegetation (isoprene, alkenes, paraffin), and (e) no lightning emissions. All model data that fall 
into the 800-400 hPa region have been used (no time sampling/averaging kernel application). 


America, an effect which is rather pronounced in the UKCA 
model. In general, we can claim that the no_VEG run pro- 
duced the most uncertain results. 

The perturbation of lightning in the noXIGHT runs 
(Fig. 8e, Fig. 9e) changes the O3-CO correlation in various 
areas in the G-PUCCINI model, like the eastern Pacific and 
the Indian Ocean, but it does not affect substantially the cor- 
relations in the UKCA model. This result initiated an exami- 
nation of the lightning parameterization in the UKCA model, 
which revealed that lightning NO x is emitted at too low tro- 
pospheric levels in the model. This may explain the low cor- 
relations shown for the UKCA model in Fig. 3 for some of 
the regions, but not everywhere, as lightning does not seem 
to be a major driver in areas like the midlatitude outflow re- 
gions (see Fig. 8e). 


5 Discussion 

Had we found an emission type with a large impact on G- 
PUCCINI correlations, but not on UKCA, and with a large 
impact on regions with high correlations (midlatitude conti- 
nental outflow regions, tropical biomass burning regions) we 
would have concluded that the treatment of particular types 
of emissions in the UKCA model (either the emissions them- 
selves or their chemical fate in the troposphere) is the reason 
for its lower correlations in important areas. Note that we 
performed an extra experiment with removed VOC emissions 
(not shown) which showed that there are not drastic changes 
in the correlations when applying such a perturbation. This 
rulled out the possibility that the different lumping of VOCs 
in the two models could make a large difference. 
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Fig. 9. Correlation coefficient of daily mean July-August 2006 middle/lower free tropospheric O3 and CO concentrations in the BASE 
UKCA run (unperturbed; a) and in 4 sensitivity runs in summer. The four runs are: (b) no anthropogenic emissions, (c) no biomass burning 
emissions, (d) no emissions from vegetation (isoprene and acetone), and (e) no lightning emissions. All model data that fall into the 800- 
400 hPa region have been used (no time sampling/averaging kernel application). 


From the above, we conclude that the discrepancies most 
likely have to do with other factors. This is supported by 
the fact that the net chemical tendencies (total O 3 chemical 
production minus total ozone chemical loss) in 80CM-00 hPa 
(Fig. 10) do not differ much in terms of sign in the crucial 
regions that we have examined (northern Atlantic, northern 
Pacific, tropical Africa). In fact, both models show net chem- 
ical O 3 destruction in large parts of the regions where we 
have found substantial positive O 3 -CO correlations. This in- 
dicates that, contrary to what is often stated, positive O 3 -CO 
correlations cannot always be used as an indication for lo- 
cal photochemical O 3 production. In areas like the north- 
ern Pacific, we expect that frequent meridional transport of 
I 0 W-CO/I 0 W-O 3 air from lower latitudes could be playing a 
central role. 


Increased STE in the UKCA model could have explained 
the lower O 3 -CO correlations, especially at midlatitudes. 
However, global STE for 2006 in the model is 659Tgyr _1 , 
which is well within reported ranges (e.g. Stevenson et al., 
2006). In fact, the global downward flux of O 3 from a surface 
(380 hPa) right above our region of interest (800-400 hPa) 
is 20-30 % lower in the UKCA model compared to the G- 
PUCCINI model. Apart from the influence of vertical mix- 
ing, we also expect that horizontal mixing can impact the cor- 
relations and their differences in the models. Additionally, 
convection can also be a source of differences. It is common 
that models forced with the same meteorology may simu- 
late vertical transport of pollution in the tropics very differ- 
ently, depending on the parameterization used (Hoyle et al., 
2010 ), while very similar (or the same) models can produce 
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Fig. 10. Net chemical tendency of middle/lower free tropospheric (800-400 hPa) O 3 in July-August 2006 for the G-PUCCINI and UKCA 
unperturbed runs (BASE). 


different results when forced with different meteorological 
analyses (Liu et al., 2010). In fact, the former study notes that 
nudged models, including UKCA, can misrepresent convec- 
tive fluxes due to misrepresentation of surface fluxes, such 
as moisture, even if the large scale circulation is well repre- 
sented. This was evidenced in too small a flux of idealised 
tracers in the UKCA model in the West African monsoon. 

Another factor that can be driving the differences between 
the models is the use of an offline photolysis scheme in 
the UKCA versus an interactive photolysis scheme in G- 
PUCCINI. As found in Voulgarakis et al. (2009a,b), the ef- 
fects of clouds on photolysis can drive sizeable changes in 
the chemistry, with impacts being stronger above and below 
clouds. Even when these effects are above/below our region 
of interest (800-400 hPa), they can still be affecting our re- 
sults via transport of signals. 

Apart from being a useful addition to future model inter- 
comparisons, the method presented here will be an invaluable 
way of evaluating the next improved version of the UKCA 
model and understanding how its recent developments (Fast- 
JX photolysis, more realistic lightning NO x ; not included 
here) impact its ability to simulate tropospheric chemistry. 

6 Conclusions 

We used TES satellite observations and two different mod- 
els to provide the first global O 3 -CO correlation analysis us- 
ing a multi-year dataset for two different seasons. A vari- 
ety of features have been found to be robust across different 
datasets and different methods of estimating the correlation, 
and also consistent with past studies that used in-situ mea- 
surements. Areas that consistently show positive correlations 
are the tropical biomass burning regions for both seasons and 
the northern midlatitude continental outflow regions in the 
summer. Among the latter, the area with the highest correla- 
tions (>0.8) is the northern Pacific, not just near east Asia, 


but even close to the west coast of North America. Areas 
of negative correlations include the high northern latitudes 
in both seasons, as well as regions like the Sahara and the 
Middle East. 

We evaluated two models’ (G-PUCCINI, UKCA) abil- 
ity to capture the O 3 -CO correlations and found that G- 
PUCCINI has a good skill, especially in summer months, but 
also has some weaknesses, especially in the winter over the 
tropics. The UKCA model, whose simulated O 3 abundances 
are not very different to these of the G-PUCCINI model, is 
able to capture some of the features (e.g. in some remote re- 
gions), but in general has correlations that are systematically 
lower than TES and G-PUCCINI over various regions. 

By performing a series of sensitivity experiments, we find 
that the only individual type of emission which is able to ex- 
plain the sign of an important correlation pattern by its own, 
is biomass burning in the tropics. The rest of the features 
of the geographical distribution remain similar when exclud- 
ing different individual types of emissions in our model ex- 
periments. However, the magnitude of the correlations do 
change, especially downwind of emission regions, where O 3 
production is affected by the perturbations that we apply. The 
analysis also reveals that emissions and chemical processing 
are most likely not the reason for the UKCA model under- 
estimates of O 3 -CO correlations. We discuss that the dis- 
crepances are more likely due to differences in vertical and 
horizontal mixing. 

There is still a lot to be examined and explained about the 
features of our correlation maps. In particular, more anal- 
ysis could be performed to understand the influence of the 
dynamics on the O 3 -CO relationship. Studying the correla- 
tions at higher altitudes would possibly enhance our under- 
standing of the role of STE and convection. It would also be 
very useful to examine the correlation patterns for other sea- 
sons, like Northern Hemisphere spring, when STE is more 
intense at midlatitudes. Upper tropospheric ozone and CO 
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data from the Aura Microwave Limb Sounder (MLS) could 
also be used for these purposes. Another interesting implica- 
tion of our approach would be to study the O 3 -CO correlation 
in a modeled future atmosphere, taking into account climate 
change and stratospheric O 3 recovery, both of which have the 
potential to influence tropospheric chemistry. 

In general, we have showed that the correlation maps pre- 
sented here are fairly robust and we believe that they are suit- 
able for use for comparisons in future model evaluations. The 
observed correlations (Fig. 3a, b) are in principle the most 
reliable, while the G-PUCCINI correlations from raw model 
output (Fig. 6 c, g) can also provide a straightforward bench- 
mark for comparison with other models (especially since raw 
model correlations have not been degraded by the observa- 
tional error). This is especially true for the Northern Hemi- 
sphere summer, when performance is even better. Over- 
all, the method presented here could be a useful addition to 
multi-model validation efforts that will be materialized in the 
future, providing insight into the models’ ability to simulate 
O 3 well for the right reasons. With further work, we may 
also be able to link discrepances in the O 3 -CO correlation 
with particular weaknesses in models, something that is not 
fully allowed by our current level of understanding. 
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